201711.02298v1 


chinaXiv 


Observation of the doubly radiative decay n’ — yyr° 


M. Ablikim!, M. N. Achasov??, S. Ahmed, X. C. Ai!, O. Albayrak, M. Albrecht*, D. J. Ambrose?, A. Amoroso?9 4.90€. 
F. F. An!, Q. An?7?5, J. Z. Bai!, O. Bakina??, R. Baldini Ferroli??^, Y. Ban?!, D. W. Bennett?®, J. V. Bennett?, 
N. Berger??, M. Bertani?°“, D. Bettoni?!^, J. M. Bian**, F. Bianchi??^99C, p, Boger?*^, I. Boyko??, R. A. Briere?, H. Cai??, 
X. Cait38, O. Cakir*^, A. Calcaterra?'^, G. F. Cao™*?, S. A. Cetin*!P, J. F. Chang'??, G. Chelkov?*^*, G. Chen!, 
H. S. Chen", J. C. Chen!, M. L. Chen! ?$, S. Chen?^, S. J. Chen??, X. Chen^?*, X. R. Chen?6, Y. B. Chen!?5, X. K. Chu?!, 
G. Cibinetto?!4, H. L. Dai'?5, J. P. Dai?^^, A. Dbeyssi'^, D. Dedovich??, Z. Y. Deng!, A. Denig??, I. Denysenko??, 
M. Destefanis?"^:907. P. De Morj?94:50€. y; Ding”, C. Dong”, J. Dong’, L. Y. Dong"? , M. Y. Dong???, Z. L. Dou?®, 
S. X. Du, P. F. Duan!, J. Z. Fan’, J. Fang! ??, S. S. Fang^?^, X. Fang? ??, Y. Fang!, R. Farinelli??4:21P, [, paya508,50C. 
F. Feldbauer??, G. Felici?" ^, C. Q. Feng^^??, E. Fioravanti?!*^, M. Fritsch??^, C. D. Fu!, Q. Gao’, X. L. Gao*7^?5, Y. Gao*9, 
Z. Gao**9*, T, Garzia?!^, K. Goetzen!?, L. Gong??, W. X. Gong!?5, W. Gradl??, M. Greco?94500 M. H. Gu?5, Y. T. Gu, 
Y. H. Guant, A. Q. Guo!, L. B. Guo?5, R. P. Guo!, Y. Guo!, Y. P. Guo??, Z. Haddadi??, A. Hafner??, S. Han??, X. Q. Hao?, 
F. A. Harris”, K. L. He!*?, F. H. Heinsius, T. Heldt, Y. K. Heng™?®*?, T. Holtmann’, Z. L. Hout, C. Hu?5, H. M. Hul*?, 
J. F. Hu994590. T. Hut3842, Y. Hut, G. S. Huang^^?5, J. S. Huang”, X. T. Huang??, X. Z. Huang??, Z. L. Huang?”, 
T. Hussain*?, W. Ikegami Andersson®!, Q. Jit, Q. P. Ji®, X. B. Jib*?, X. L. Ji?5, L. W. Jiang??, X. S. Jiang! ?94?, 

X. Y. Jiang’, J. B. Jiao?*, Z. Jiao!*, D. P. Jin^?94, S. Jint*?, T. Johansson®', A. Julin**, N. Kalantar-Nayestanaki??, 
X. L. Kang!, X. S. Kang®°, M. Kavatsyuk??, B. C. Keř, P. Kiese??, R. Kliemt!?, B. Kloss??, O. B. Kolcu** P^, B. Kopf*, 
M. Kornicer?, A. Kupsc?', W. Kiihn?*, J. S. Lange?*, M. Lara?^, P. Larin'^, H. Leithoff??, C. Leng?’ , C. Li?!, 
Cheng Li*^^?5, D. M. Li?*, F. Lit?8, F. Y. Li?! , G. Li!, H. B. Li^??, H. J. Lit, J. C. Lit, Jin Li?, K. Li?, K. Lit’, Lei Li?, 
P. R. Li?^*, Q. Y. Li?, T. Li?, W. D. Li^?, W. G. Li!, X. L. Li?, X. N. Li. X. Q. Li?^, v. B. Li?, Z. B. Li??, 

H. Liang^^?5, Y. F. Liang”, Y. T. Liang?*, G. R. Liao", D. X. Lin!^, B. Liu?^^, B. J. Liu!, C. X. Liu!, D. Liu*^?5, 

F. H. Liu??, Fang Liu’, Feng Liu?, H. B. Liu?, HuanHuan Liu!, HuiHui Liu!9$, H. M. Liu’, J. Liu!, J. B. Liu?^?5, 

J. P. Liu??, J. Y. Liu’, K. Liu?^?, K. Y. Liu?’, L. D. Liu?!, P. L. Liu*?5, Q. Liu, S. B. Liu4^?5, X. Liu?6, Y. B. Liu®?, 
Y. Y. Liu??, Z. A. Liu, Zhiqing Liu??, H. Loehner??, X. C. Lou!??:2, H. J. Lut’, J. G. Lu, Y. Lut, Y. P. Lut”, 
C. L. Luo??, M. X. Luo??, T. Luo*?, X. L. Luo^?5, X. R. Lyu?, F. C. Ma”, H. L. Mat, L. L. Ma®?, M. M. Mat, Q. M. Ma! 
T. Mat, X. N. Ma”, X. Y. Mal, y. M. Ma’, F. E. Maas!*, M. Maggiora?'^:??7. Q. A. Malik’, Y. J. Mao?! , Z. P. Maot, 
S. Marcello??^:5997. J. G. Messchendorp??, G. Mezzadri?!?, J. Min!?5, T. J. Mint, R. E. Mitchell?, X. H. Mo!?5:22, 

Y. J. Mo®, C. Morales Morales!^, N. Yu. Muchnoi??, H. Muramatsu**, P. Musiol*, Y. Nefedov??, F. Nerling!?, 

I. B. Nikolaev?:?, Z. Ning'?5, S. Nisar?, S. L. Niu^?5, X. Y. Niu!, S. L. Olsen??, Q. Ouyang! ??:2?, S. Pacetti?"P, Y, Pan?*^?5, 
M. Papenbrock?!, P. Patteri??^, M. Pelizaeus^, H. P. Peng^^35, K. Peters!?9?, J. Pettersson?! , J. L. Ping”, R. G. Ping^??, 
R. Poling**, V. Prasad!, H. R. Qi?, M. Qi”, S. Qian!?$, C. F. Qiao?, L. Q. Qin??, N. Qin??, X. S. Qin!, Z. H. Qin ?5, 

J. F. Qiu!, K. H. Rashid ???, C. F. Redmer??, M. Ripka”, G. Rong!*?, Ch. Rosner!^, X. D. Ruan??, A. Sarantsev?**, 

M. Savri&?!P, C. Schnier*, K. Schoenning?!, W. Shan?!, M. Shao^^?5, C. P. Shen?, P. X. Shen’, X. Y. Shen ??, 

H. Y. Sheng’, W. M. Song!, X. Y. Song’, S. Sosio?94:59C. S. Spataro? ^97. G, X. Sunt, J. F. Sun?®, S. S. Sunt*?, 
. H. Sunt, Y. J. Sun**?5, Y. Z. Sun! , Z. J. Sun'??, Z. T. Sun??, C. J. Tang®®, X. Tang!, I. Tapan*'^, E. H. Thorndike?, 
Tiemens??, I. Uman* P, G. S. Varner?, B. Wang??, B. L. Wang*”, D. Wang?! , D. Y. Wang?!, K. Wang'??, L. L. Wang", 
Wang!, M. Wang??, P. Wang!, P. L. Wang!, W. Wang'??, W. P. Wang*”°8, X. F. Wang^^, Y. Wang?', Y. D. Wang", 
. Wang! ?52?, Y. Q. Wang”, Z. Wang, Z. G. Wang'?5, Z. H. Wang '^?5, Z. Y. Wang!, Zongyuan Wang!, T. Weber??, 
. Wei!!, P. Weidenkaff??, S. P. Went, U. Wiedner*, M. Wolke?!, L. H. Wut, L. J. Wut, Z. Wu, L. Xia7^?5, L. G. Xia“, 


, 


S. 
F 
H 
Xia!®, D. Xiao!, H. Xiao*?, Z. J. Xiao??, Y. G. Xie! ?5, Y. H. Xie, Q. L. Xiu!?5, G. F. Xut, J. J. Xut, L. Xut, Q. J. Xu”, 
N. Xu? , X. P. Xu?*, L. Yan?94900. w. B. Yan755, W. C. Yan!7?*, Y. H. Yan!5, H. J. Yang?^^, H. X. Yang!, L. Yang??, 
Y. X. Yang!!, M. Ye!?5, M. H. Ye", J. H. Yin!, Z. Y. You??, B. X. Yut3®4?, C. X. Yu?^, J. S. Yu?5, C. Z. Yuan!??, 
Y. Yuan!, A. Yuncu*!?:*, A. A. Zafar?, Y. Zeng!®, Z. Zeng^'?5, B. X. Zhang!, B. Y. Zhang'?5, C. C. Zhang", 
D. H. Zhang! , H. H. Zhang??, H. Y. Zhang'??, J. Zhang! , J. J. Zhang! , J. L. Zhang! , J. Q. Zhang!, J. W. Zhang! ?9:4?, 
J. Y. Zhang’, J. Z. Zhang", K. Zhang!, L. Zhang!, S. Q. Zhang’, X. Y. Zhang??, Y. Zhang!, Y. H. Zhang! ?5, 
Y. N. Zhang? , Y. T. Zhang?^^?*, Yu Zhang? , Z. H. Zhang?, Z. P. Zhang”, Z Y. Zhang??, G. Zhao!, J. W. Zhao!?5, 
J. Y. Zhao!, J. Z. Zhao'??, Lei Zhao*”**, Ling Zhao!, M. G. Zhao??, Q. Zhao!, Q. W. Zhao!, S. J. Zhao?*, T. C. Zhao!, 
Y. B. Zhao! ?$, Z. G. Zhao*^??, A. Zhemchugov??^, B. Zheng^*' ^, J. P. Zheng! ??, W. J. Zheng??, Y. H. Zheng*?, 
B. Zhong??, L. Zhou??, X. Zhou??, X. K. Zhou**?5, X. R. Zhou*^??, X. Y. Zhou!, K. Zhu!, K. J. Zhu! ?9/4?, S. Zhu!, 
S. H. Zhu*®, X. L. Zhu??, Y. C. Zhu^^?5, Y. S. Zhut4?, Z. A. Zhu^?, J. Zhuang'?*, L. Zotti?9^99C. B, S. Zou!, J. H. Zou! 


(BESHI Collaboration) 


1 Institute of High Energy Physics, Beijing 100049, People's Republic of China 
? Beihang University, Beijing 100191, People's Republic of China 
? Beijing Institute of Petrochemical Technology, Beijing 1026177, People's Republic of China 
^ Bochum Ruhr-University, D-44780 Bochum, Germany 
5 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
$ Central China Normal University, Wuhan 430079, People’s Republic of China 
T China Center of Advanced Science and Technology, Beijing 100190, People’s Republic of China 
8 COMSATS Institute of Information Technology, Lahore, Defence Road, Off Raiwind Road, 54000 Lahore, Pakistan 
? G.I. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia 
19 GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany 


x 
M. 
in 
Y. 
D. 
y. 
Q. 


201711.02298v1 


chinaXiv 


11 Guangxi Normal University, Guilin 541004, People's Republic of China 
1 Guangxi University, Nanning 530004, People’s Republic of China 
13 Hangzhou Normal University, Hangzhou 310036, People’s Republic of China 
14 Helmholtz Institute Mainz, Johann-Joachim-Becher-Weg 45, D-55099 Mainz, Germany 
15 Henan Normal University, Xinziang 453007, People’s Republic of China 
16 Henan University of Science and Technology, Luoyang 471003, People's Republic of China 
17 Huangshan College, Huangshan 245000, People’s Republic of China 
18 Hunan University, Changsha 410082, People’s Republic of China 
19? Indiana University, Bloomington, Indiana 47405, USA 
20 (A)INFN Laboratori Nazionali di Frascati, I-00044, Frascati, 
Italy; (B)INFN and University of Perugia, I-06100, Perugia, Italy 
21 (A)INFN Sezione di Ferrara, I-44122, Ferrara, Italy; (B)University of Ferrara, I-44122, Ferrara, Italy 
22? Johannes Gutenberg University of Mainz, Johann-Joachim-Becher- Weg 45, D-55099 Mainz, Germany 
23 Joint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia 
24 Justus-Liebig-Universitaet Giessen, II. Physikalisches Institut, Heinrich-Buff-Ring 16, D-35392 Giessen, Germany 
25 KVI-CART, University of Groningen, NL-9747 AA Groningen, The Netherlands 
26 Lanzhou University, Lanzhou 730000, People’s Republic of China 
27 Liaoning University, Shenyang 110036, People’s Republic of China 
?5 Nanjing Normal University, Nanjing 210023, People’s Republic of China 
?? Nanjing University, Nanjing 210093, People’s Republic of China 
30 Nankai University, Tianjin 300071, People’s Republic of China 
?! Peking University, Beijing 100871, People’s Republic of China 
32 Seoul National University, Seoul, 151-747 Korea 
33 Shandong University, Jinan 250100, People’s Republic of China 
34 Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China 
35 Shanxi University, Taiyuan 030006, People’s Republic of China 
36 Sichuan University, Chengdu 610064, People’s Republic of China 
57 Soochow University, Suzhou 215006, People’s Republic of China 
38 State Key Laboratory of Particle Detection and Electronics, Beijing 100049, Hefei 230026, People’s Republic of China 
39 Sun Yat-Sen University, Guangzhou 510275, People’s Republic of China 
^9 Tsinghua University, Beijing 100084, People’s Republic of China 
^! (A)Ankara University, 06100 Tandogan, Ankara, Turkey; (B)Istanbul Bilgi University, 34060 Eyup, Istanbul, Turkey; 
(C)Uludag University, 16059 Bursa, Turkey; (D)Near East University, Nicosia, North Cyprus, Mersin 10, Turkey 
? University of Chinese Academy of Sciences, Beijing 100049, People’s Republic of China 
^5 University of Hawaii, Honolulu, Hawaii 96822, USA 
44 University of Minnesota, Minneapolis, Minnesota 55455, USA 
^ University of Rochester, Rochester, New York 14627, USA 
^6 University of Science and Technology Liaoning, Anshan 114051, People’s Republic of China 
47 University of Science and Technology of China, Hefei 230026, People’s Republic of China 
5 University of South China, Hengyang 421001, People’s Republic of China 
49 University of the Punjab, Lahore-54590, Pakistan 
50 (A)University of Turin, I-10125, Turin, Italy; (B)University of Eastern 
Piedmont, 1-15121, Alessandria, Italy; (C)INFN, 1-10125, Turin, Italy 
9! Uppsala University, Box 516, SE-75120 Uppsala, Sweden 
52 Wuhan University, Wuhan 430072, People’s Republic of China 
53 Zhejiang University, Hangzhou 310027, People’s Republic of China 
54 Zhengzhou University, Zhengzhou 450001, People’s Republic of China 


^ Also at Bogazici University, 34342 Istanbul, Turkey 
^ Also at the Moscow Institute of Physics and Technology, Moscow 141700, Russia 
^ Also at the Functional Electronics Laboratory, Tomsk State University, Tomsk, 634050, Russia 
7 Also at the Novosibirsk State University, Novosibirsk, 630090, Russia 
^ Also at the NRC "Kurchatov Institute", PNPI, 188300, Gatchina, Russia 
^ Also at Istanbul Arel University, 34295 Istanbul, Turkey 
9? Also at Goethe University Frankfurt, 60323 Frankfurt am Main, Germany 


^ Also at Key Laboratory for Particle Physics, Astrophysics and Cosmology, Ministry of Education; Shanghai Key Laboratory 
for Particle Physics and Cosmology; Institute of Nuclear and Particle Physics, Shanghai 200240, People's Republic of China 


* Government College Women University, Sialkot - 51310. Punjab, Pakistan. 


Based on a sample of 1.31 billion J/w events collected with the BESIII detector, we report the 
study of the doubly radiative decay n’ — yyr° for the first time, where the 7’ meson is produced via 
the J/4» — yn! decay. The branching fraction of n! — yyr° inclusive decay is measured to be B(r' > 
VT" Jma. = (3.20 + 0.07(stat) + 0.23(sys)) x 107°, while the branching fractions of the dominant 
process n’ — yw and the nonresonant component are determined to be B(y' — yw) x B(w — yr?) = 
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(23.7 + 1.4(stat) + 1.8(sys)) x 1074 and B(n’ + yyr°)nwr = (6.16 + 0.64(stat) + 0.67(sys)) x 1074, 
respectively. In addition, the M32, -dependent partial widths of the inclusive decay are also presented. 


PACS numbers: 13.40.Gp, 13.40.Hq, 13.20.Jf, 14.40.Be 


I. INTRODUCTION 


The 7 meson provides a unique stage for understand- 
ing the distinct symmetry-breaking mechanisms present 
in low-energy quantum chromodynamics (QCD) 
and its decays play an important role in exploring the 
effective theory of QCD at low energy |6]. Recently, the 
doubly radiative decay 7! — yyr° was studied in the 
frameworks of the linear c model (LoM) and the vec- 
tor meson dominance (VMD) model [4 [8]. It has been 
demonstrated that the contributions from the VMD are 
dominant. Experimentally, only an upper limit of the 
nonresonant branching fraction of B(n! > yy7°)nr < 
8x 1074 at the 90% confidence level has been determined 
by the GAMS-2000 experiment id. 


In this article, we report the first measurement of the 
branching fraction of the inclusive 7’ — yyr° decay and 
the determination of the M. dependent partial widths, 
where M, is the invariant mass of the two radiative pho- 
tons. The inclusive decay is defined as the 7’ decay into 
the final state yy” including all possible intermediate 
contributions from the p and w mesons below the n’ mass 
threshold and the nonresonant contribution from the ex- 
cited vector meson above the 7 mass threshold. Since 
the contribution from mesons above the 7’ threshold ac- 
tually derives from the low-mass tail and looks like a con- 
tact term, we call this contribution ‘nonresonant’. The 
branching fraction for the nonresonant n! — yyr° decay 
is obtained from a fit to the yr° invariant mass distribu- 
tion by excluding the coherent contributions from the p 
and w intermediate states. The measurement of the M2, 
dependent partial widths will provide direct inputs to the 
theoretical calculations on the transition form factors of 
n! — yyr° and improve the theoretical understanding of 
the 7 decay mechanisms. 


Il. EXPERIMENTAL DETAILS 


The source of 7’ mesons is the radiative J/W — yn’ 
decay in a sample of 1.31 x 10° J/w events Ed roi col- 
lected by the BESIII detector. Details on the features 
and capabilities of the BESIII detector can be found in 
Ref. r 


'The response of the BESIII detector is modeled with 
a Monte Carlo (MC) simulation based on GEANT4 [13]. 
The program EVTGEN is used to generate a J/y — 
yn’ MC sample with an angular distribution of 1+cos? 04, 
where 0, is the angle of the radiative photon relative to 
the positron beam direction in the J/4» rest frame. The 
decays n! — yw(p), w(p) + "yn? are generated using the 


helicity amplitude formalism. For the nonresonant n’ + 
"yy? decay, the VMD model TIE) is used to generate the 
MC sample with p(1450) or w(1650) exchange. Inclusive 
J/w decays are generated with KKMC generator; the 
known J/w decay modes are generated by EVTGEN 
with branching fractions setting at Particle Data Group 
(PDG) world average values |16]; the remaining unknown 
decays are generated with LUNDCHARM [1]. 


IIl. EVENT SELECTION AND BACKGROUND 
ESTIMATION 


Electromagnetic showers are reconstructed from clus- 
ters of energy deposits in the electromagnetic calorime- 
ter (EMC). The energy deposited in nearby time-of-light 
(TOF) counters is included to improve the reconstruction 
efficiency and energy resolution. The photon candidate 
showers must have a minimum energy of 25 MeV in the 
barrel region (|cos0| < 0.80) or 50 MeV in the end cap 
region (0.86 « |cos0| « 0.92). Showers in the region 
between the barrel and the end caps are poorly mea- 
sured and excluded from the analysis. In this analysis, 
only the events without charged particles are subjected 
to further analysis. The average event vertex of each run 
is assumed as the origin for the selected candidates. To 
select J/1 — yn, n! — yyr° (19 — yy) signal events, 
only the events with exactly five photon candidates are 
selected. 


To improve resolution and reduce background, a 
five-constraint kinematic (5C) fit imposing energy- 
momentum conservation and a 7° mass constraint is per- 
formed to the yyy7° hypothesis, where the 1? candidate 
is reconstructed with a pair of photons. For events with 
more than one 7° candidate, the combination with the 
smallest y2¢ is selected. Only events with xe < 30 are 
retained. The x26 distribution is shown in Fig. [I] with 
events in the r/ signal region of |M;,,o — My | < 25 MeV 
(M is the 7’ nominal mass from PDG [18]). In order to 
suppress the multi-7° backgrounds and remove the mis- 
combined 7° candidates, an event is vetoed if any two of 
five selected photons (except for the combination for the 
n° candidate) satisfies |My. — Mzo| < 18 MeV/c?, where 
Mpo is the 7° nominal mass. After the application of the 
above requirements, the most energetic photon is taken 
as the primary photon from the J/w decay, and the re- 
maining two photons and the 7° are used to reconstruct 
the ņ candidates. Figure 2] shows the yyr° invariant 
mass spectrum. 


Detailed MC studies indicate that no peaking back- 
ground remains after all the selection criteria. The 
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FIG. 1: Distribution of the x26 of the 5C kinematic fit for 
the inclusive 7’ decay. Dots with error bars are data; the 
heavy (black) solid-curve is the sum of signal and expected 
backgrounds from MC simulations; the light (red) solid-curves 
is signal components which are normalized to the fitted yields; 
the (green) dotted-curve is the class I background; and the 
(pink) dot-dashed-curve is the class II background. 
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FIG. 2: Results of the fit to Myro for the selected inclusive 


n’ > yyr° signal events. The (black) dots with error bars are 
the data. 


sources of backgrounds are divided into two classes. 
Background events of class I are from J/v — yn! with 
7 decaying into final states other than the signal final 
states. These background events accumulate near the 
lower side of the 7’ signal region and are mainly from 
n — wr n (n > yy), y! 2 37? and 1! — yy, as shown 
as the (green) dotted curve in Fig. [2] Background events 
in class II are mainly from J/~w decays to final states with- 
out 7’, such as J/ => yr°r® and J/y — wn (w > yr, 
1] — yy) decays, which contribute a smooth distribution 
under the 7’ signal region as displayed as the (pink) dot- 
dashed curve in Fig. [2] 


IV. SIGNAL YIELDS AND BRANCHING 
FRACTIONS 


A fit to the yyr° invariant mass distribution is per- 
formed to determine the inclusive n — yyr° signal yield. 
The probability density function (PDF) for the signal 
component is represented by the signal MC shape, which 
is obtained from the signal MC sample generated with an 
incoherent mixture of p, w and the nonresonant compo- 
nents according to the fractions obtained in this analysis. 
Both the shape and the yield for the class I background 
are fixed to the MC simulations and their expected inten- 
sities. The shape for the class II background is described 
by a third-order Chebychev polynomial, and the corre- 
sponding yield and PDF parameters are left free in the 
fit to data. The fit range is 0.70—1.10 GeV/c?. Figure P] 
shows the results of the fit. The fit quality assessed with 
the binned distribution is y?/n.d.f = 108/95 = 1.14. The 
signal yield and the MC-determined signal efficiency for 
the inclusive n’ decay are summarized in Table [I] 


In this analysis, the partial widths can be obtained 
by studying the efficiency-corrected signal yields for each 
given M2, bin i for the inclusive n — yyr° decay. The 
resolution in M2, is found to be about 5x 10? (MeV /c?)? 
from the MC simulation, which is much smaller than 
1.0 x 10* (MeV /c?)?, a statistically reasonable bin width, 
and hence no unfolding is necessary. The 7’ signal yield in 
each M2, bin is obtained by performing bin-by-bin fits to 
the yz? invariant mass distributions using the fit proce- 
dure described above. Thus the background-subtracted, 
efficiency-corrected signal yield can be used to obtain the 
partial width for each given M, R interval, where the PDG 
value is used for the total width of the 7’ meson ie. The 
results for dP(n' — yy7?)/d M2, in each MZ, interval are 
listed in Table [IT] and depicted in Fig. [3] where the con- 
tributions from each component obtained from the MC 
simulations are normalized with the yields by fitting to 
M,,0 as displayed in Fig. 


Assuming that the inclusive decay n! > yyr° can be 
attributed to the vector mesons p and w and the nonres- 
onant contribution, we apply a fit to the yr° invariant 
mass to determine the branching fraction for the non- 
resonant n! — yyr” decay using the 7’ signal events 
with |My 0 — m,,| < 25 MeV/c?. In the fit, the p- 
w interference is considered, but possible interference 
between the w (p) and the nonresonant process is ne- 
glected. To validate our fit, we also determine the prod- 
uct branching fraction for the decay chain 7’ — yw, 
w — yn”. Figure[l]shows the Mo distribution. Since 
the doubly radiative photons are indistinguishable, two 
entries are filled into the histogram for each event. For 
the PDF of the coherent w and p produced in 7! — 
yyn°, we use [e(My_0) x B x E da x [BW (Myro) + 
ae? BW, (M. o)? x B5, xB2(]9G(0, v), where e(M,,0) 
is the detection efficiency determined by the MC simula- 
tions; E, ,/(,/5) iS the energy of the transition photon in 


201711.02298v1 


chinaXiv 


TABLE I: Observed 7 signal yields (N n^) and detection efficiencies (c) for inclusive n > yyr°, n > yw(w — ^yn?), and the 
nonresonant n — yyr° decays. The measured branching fractions? in this work, comparison of values from the PDG and 
theoretical predictions are listed. The first errors are statistical and the second ones are systematic. 


n! > yyr° (Inclusive) n! + yw,w — yr? n! — yz? (Nonresonant) 


7 


N” 3435 + 76 + 244 
€ 16.1% 
B (1074) 32.0 + 0.7 + 2.3 
Brena (1074) 


Predictions (1074) 57 [7],65 [8] 


2340 + 141 + 180 655 + 68 + 71 
14.8% 15.9% 
23.7 1.4 1.8? 6.16 + 0.64 + 0.67 

21.7 x 1.3" «8 


a The product branching fraction B(n! — yw) - B(w — yr°). © The product branching fraction B(n! > yw) - B(w — y7?) from PDG Ed. 
© The product branching fraction B(n! — yp?) - B(p9 — y7?) is determined to be (1.92 + 0.16(stat)) x 1074 using the fitted yield in 
Fig. Ø] which is in agreement with the PDG value of (1.75 + 0.23) x 1074 [16]. 


TABLE II: Results for dr (n! — yyr°)/dM?3, (in units of keV/(GeV /c?)?) for thirteen intervals of M?,. The first uncertainties 


are statistical and the second systematic. 


M2, ((GeV/c?)?) 
aD (n' 2 yyn?)/M2 
M3, ((GeV/c?)?) 
dU(n' — yyn?)/M$ 

M32., ((GeV/c?)?) 
dU(' 2 yyn?)/M$ 


- 0.18 d 


the rest frame of 7 (w/p); BW..(M,,0) is a relativistic 
Breit-Wigner (BW) function, and BW, (M. zo) is a rel- 
ativistic BW function with mass-dependent width is. 
The masses and widths of the p and w meson are fixed to 
their PDG values [16]. B2, /p) 18 the Blatt-Weisskopf 
centrifugal barrier factor for the m'(w/p) decay vertex 
with radius R = 0.75 fm [19, 2d], and B2,,,,;,; is used to 
damp the divergent tail due to the factor ET whe)" The 


Gaussian function G(0, v) is used to parameterize the de- 
tector resolution. The combinatorial background is pro- 
duced by the combination of the 7° and the photon from 
the 7’ meson, and its PDF is described with a fixed shape 
from the MC simulation. The ratio of yields between the 
combinatorial backgrounds and the coherent sum of p-w 
signals is fixed from the MC simulations. The shape of 
the nonresonant signal 7’ — ^a? is determined from the 
MC simulation, and its yield is determined in the fit. The 
background from the class I as discussed above is fixed to 
the shape and yield of the MC simulation. Finally, the 
shape from the class IT background is obtained from the 
n’ mass sidebands (738—788 and 1008—1058 MeV/c?), 
and its normalization is fixed in the fit. The Myro mass 
range used in the fit is 0.20—0.92 GeV /c?. In the fit, the 
interference phase 0 between the p- and w-components is 
allowed. Due to the low statistics of the p meson con- 
tribution, we fix the ratio a of p and w intensities to 
the value for the ratio of B(n’ — yp) - B(p — yT?) and 
B(n! — qw) - B(w — yr°) from the PDG [16]. Figure H] 


[0.04, 0.06] [0.06, 0.09] (0.09, 0.12] 
2.60 + 0.15 + 0.18 1.87 + 0.12 + 0.14 | 1.76 + 0.11 + 0.13 
[0.20, 0.25] 0.25, 0.28] 
1.97 + 0.10 + 0.14 |2.00 + 0.17 + 0.15 |1.07 + 0.20 + 0.08 
] 


5 
[0 4 


0. 


shows the results. The yields for the vector mesons p, w 
and their interference are determined to be (183 + 15), 
(2340 + 141), and (174 + 92), respectively. The sig- 
nal yields and efficiencies as well as the corresponding 
branching fractions for the n! — yw(w — ^x?) and non- 
resonant decays are summarized in Table [I] 


V. SYSTEMATIC UNCERTAINTIES 


'The systematic uncertainties on the branching fraction 
measurements are summarized in Table 'The uncer- 
tainty due to the photon reconstruction is determined to 
be 196 per photon as described in Ref. [1]. 'The uncer- 
tainties associated with the other selection criteria, kine- 
matic fit with x2, < 30, the number of photons equal to 
5 and 7° veto (|My — Mzo| > 18 MeV/c?) are studied 
with the control sample J/1» — yn’, n — yw, w — yr? 
decay, respectively. The systematic error in each of the 
applied selection criteria is numerically estimated from 
the ratio of the number of events with and without the 
corresponding requirement. The corresponding resulting 
efficiency differences between data and MC (2.796, 0.596, 
and 1.996 , respectively) are taken to be representative of 
the corresponding systematic uncertainties. 


In the fit for the inclusive 7’ decay, the signal shape 
is fixed to the MC simulation. The uncertainty due to 
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FIG. 3: Partial width (in keV) versus M2, for the inclusive 
n! — yyn° decay. The error includes the statistic and sys- 
tematic uncertainties. The (blue) histogram is the sum of 
an incoherent mixture of p-w and the nonresonant compo- 
nents from MC simulations; the (back) dotted-curves is w- 
contribution; the (red) dot-dashed-curve is the p-contribution; 
and the (green) dashed-curve is the nonresonant contribution. 
All the components are normalized using the yields obtained 
in Fig. 
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FIG. 4: Distribution of the invariant mass M,,o and fit 
results in the 7’ mass region. The points with error bars are 
data; the (black) dotted-curve is from the w-contribution; the 
(red) long dashed-curve is from the p-contribution; the (blue) 
short dashed-curve is the contribution of p-w interference; the 
(green) long dashed curve is the nonresonance; the (pink) 
histogram is from the class II background; the (black) short 
dot-dashed curve is the combinatorial backgrounds of 7’ > 
yw, yp. The (blue) solid line shows the total fit function. 


the signal shape is considered by convolving a Gaussian 
function to account for the difference in the mass resolu- 
tion between data and MC simulation. In the fit to the 
yn? distribution, alternative fits with the mass resolution 
left free in the fit and the radius R in the barrier factor 
changed from 0.75 fm to 0.35 fm are performed, and the 


changes of the signal yields are taken as the uncertainty 
due to the signal shape. 


In the fit to the M. o distribution, the signal shape is 
described with an incoherent sum of contributions from 
processes involving p and w and nonresonant processes 
obtained from MC simulation, where the nonresonant 
process is modeled with the VMD model. A fit with an 
alternative signal model for the different components, i.e. 
a coherent sum for the p-, w-components and a uniform 
angular distribution in phase space (PHSP) for the non- 
resonant process, is performed. The resultant changes in 
the branching fractions are taken as the uncertainty re- 
lated to the signal model. An alternate fit to the Myro 
distribution is performed, where the PDF of the nonreso- 
nant decay is extracted from the PHSP MC sample. The 
changes in the measured branching fractions are consid- 
ered to be the uncertainty arising from the signal model. 


In the fit to the M,,0 distribution, the uncertainty due 
to the fixed relative p intensity is evaluated by changing 
its expectation by one standard deviation. An alterna- 
tive fit in which the ratio of yields between combinato- 
rial backgrounds and the coherent sum of p — w signals is 
changed by one standard deviation from the MC simula- 
tion is performed, and the change observed in the signal 
yield is assigned as the uncertainty. A series of fits us- 
ing different fit ranges is performed and the maximum 
change of the branching fraction is taken as a systematic 
uncertainty. 


The uncertainty due to the class I background is esti- 
mated by varying the numbers of expected background 
events by one standard deviation according to the errors 
on the branching fraction values in PDG . The un- 
certainty due to the class II background is evaluated by 
changing the order of the Chebychev polynomial from 3 
to 4 for the fit to the n’ inclusive decay, and varying the 
ranges of 7 sidebands for the fit to the yx° invariant 
mass distribution, respectively. 


The number of J/1» events is N;,,, = (1310.6 + 10.5) x 
106 [10, [11], corresponding to an uncertainty of 0.8%. 
The branching fractions for the J/1» — yn! and 7° > yy 
decays are taken from the PDG , and the correspond- 
ing uncertainties are taken as a systematic uncertainty. 
The total systematic errors are 7.196, 7.7%, 10.8% for 
the inclusive decay, w contribution and nonresonant de- 
cay, respectively, as summarized in Table [III] 


VI. SUMMARY 


In summary, with a sample of 1.31 x 10? J/1/ events col- 
lected with the BESIII detector, the doubly radiative de- 
cay n! + yyr° has been studied. The branching fraction 
of the inclusive decay is measured for the first time to be 
B(n! — yyn?)ma. = (3.20-£0.07 (stat) +0.23(sys)) x 1073. 
The M2, dependent partial decay widths are also de- 
termined. In addition, the branching fraction for the 
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TABLE III: Summary of relative systematic uncertainties (%) 
for the branching fraction measurements. Here ninc., 1j; and 
"Nn, represent the inclusive n — yyr°, n > yw(w > y7”) 
and nonresonant decays, respectively. 


Minal. Nw NNR. 
Photon detection 
5C kinematic fit 
Number of photons 
n° veto 
Signal shape 
Signal model 
p relative intensity 
Combinatorial backgrounds 
Fit range 
Class I background 
Class II background 
Cited branching fractions 
Number of J/1» events 


Total systematic error 


nonresonant decay is determined to be B(r' > yyr°)nr 
= (6.16 + 0.64(stat) + 0.67(sys)) x 1074, which agrees 
with the upper limit measured by the GAMS-2000 ex- 
periment As a validation of the fit, the product 
branching fraction with the omega intermediate state in- 
volved is obtained to be B(r' > qw): B(w — yr°) = 
(2.37 +0.14(stat) + 0.18(sys)) x 1073, which is consistent 
with the PDG value [1g]. These results are useful to test 
QCD calculations on the transition form factor, and pro- 
vide valuable inputs to the theoretical understanding of 
the light meson decay mechanisms. 
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